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The reaction ttN-^ujN revisited: the a;— A^ scattering length 
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We reinvestigate the experimental data on tj-meson production in pion-nucleon collisions lil, 
H, m] based on both an analytical approach and a Monte Carlo simulation. Our analysis allows 
us to study the kinematical peculiarities of the 7r~p— »wn reaction at energies close to threshold. 
Based on the hypothesis, that the formalism derived in ref. |l| was applied improperly, the 
unusually strong energy dependence claimed for the near-threshold cross section is identified 
as a purely kinematical effect. Based on this hypothesis, we deduce an effective n'p^um cross 
section that is larger by an order of magnitude compared to the one commonly used. In addition 
we extract a lower limit for the imaginary part QaujN of the uiN scattering length. We find 
a value of Q'atjjv=0.24±0.05 fm, which is significantly larger than that presently found in the 
literature. 
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I. INTRODUCTION 

In Refs. 0, 0, H the reaction n^p — > con was measured 
at energies close to the nominal to production thresh- 
old. In Ref. B it was claimed that these experiments 
were interpreted incorrectly. In this work we reinvesti- 
gate the experiments using an event-by-event simulation 
as well as an analytical calculation. The purpose of this 
manuscript is two fold: Besides studying, if the reaction 
kinematics was treated properly in Refs. ||l|, ||, ||], we 
also take the oportunity to discuss in detail the method 
proposed in Ref. M for the near threshold production of 
narrow resonances, for it might prove useful! for upcom- 
ing experiments at modern accelerators. 

The reaction ttN^ujN, especially at energies close to 
the reaction threshold, attracts great interest for several 
reasons. 

1) There is a large number of baryonic resonances that 
have been predicted theoretically but not yet observed 
experimentally [|6| . Since almost all our knowledge about 
resonances is deduced from elastic ttN scattering, to re- 
ally resolve the issue of the missing resonances other 
channels must be studied. There are predictions of a 
number of baryonic resonances with masses from 1.8 to 
2 GeV which might couple to the ujN channel R|. The 
■nN—^ioN data is an essential ingredient in any kind of 
analysis aimed at identifying baryonic resonances possi- 
bly coupled to the w-meson. A reliable extraction of the 
resonance properties implies the analysis of pion and pho- 
ton induced data simultaneously for as many final states 
as possible ||, 10|. 

In the experiments described in Refs. |I[ || || the 
squared matrix element \M\^ was extracted from a mea- 
surement of the reaction vr^p — > ujn. The authors claim 
to have measured a very strong near-threshold suppres- 
sion, which might stem from the production of a baryonic 



resonance in a p-wave. The measured differential cross 
sections, however, are isotropic. Up to now a suppression 
of the matrix element as reported in p, 0, 0] cannot be 
understood theoretically |^. In any case, the 7r~p— >wn 
cross section is very important for resolving the existence 
of resonances not yet observed. 

2) The TT~p^ujn data 0] are intensively used in 
the evaluation of the in-medium properties the of tu- 
meson O, Ol O, 113. The in-medium w-meson mass 
and width are commonly considered either for low finite 
uj momenta Puj or at Pu/=0- Within the tp approximation 
an additional in-medium collisional width AT^^ of the uj- 
meson is given by the imaginary part of the forward ujN 
scattering amplitude 3/a;(0) as |ll, HQ, [l7[. 
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where ttin and rriaj are the free nucleon and w-meson 
masses, respectively, while pg is the baryon density. 
^fui{0) can be evaluated using the optical theorem from 
the total ojN cross section a^^ as 



^UO)^^<yl%, 



(2) 



with p^ denoting the w-meson momentum. 

The partial luN^ttN cross section can be obtained 
from the ttN^ujN data using detailed balance. The total 
wN cross section consists of several different final chan- 
nels, but the luN^'TtN channel is believed to be the one 
numerically most important. In any case, knowledge of 
this inelastic channel provides at least a lower bound for 
the imaginary part of the forward scattering amplitude 
3/ij(0). The above procedure was extensively applied in 
the hterature 0, |l2) . 

3) The feasibility of producing nuclear bound states of 
w-mesons |l8, nol can be estimated reliably in terms of 



the effective uiN scattering length, which again is dom- 
inantly given by the 'K~p^ujn data near the reaction 
threshold ||T^. Thus the data given in Ref. g] directly 
lead to an estimate of the possible existence of w-mesonic 
nuclei. 

In addition, it should be stressed that the near- 
threshold production of the w-mesons in the pd—^Lu^He 
reaction po| indicates that the squared matrix element 
\A4\^ strongly decreases starting from the center-of-mass 
w-meson momentum q^:^120 MeV/c to the reaction 
threshold where q^j—O- Within the indicated q^^ range 
the reduction of \Ai\'^ accounts for a factor of almost 4. 
As we discussed above, exactly the same reduction within 
the same range of the w-meson momenta quj<120 MeV/c 
was claimed for the 7r~p— s-wn measurements |3|. 

It should be stressed, however, that the squared matrix 
element |A^P evaluated from the data on other heavy 
mesons as 77' |2l| and p2| production in ttN reactions 
did not show such a strong near-threshold suppression. 
This difference might be ascribed to the fact that the 77' 
and 0-meson widths are significantly smaller than that 
of the uj. 

The above motivation still does not cover all subjects 
in which 7r~p— >u;n data were used. In a recent publi- 
cation it was argued that the experiments as described, 
e.g., in Ref. Q] were misinterpreted [^. The authors 
proposed a correction factor that completely removed 
the above-mentioned strong energy dependence from the 
cross section. If confirmed, this factor would answer why 
all microscopic calculations performed so far |g^ failed to 
reproduce the energy dependence of the data correctly. 
The main conclusion of this work will be to argue, that, 
although the formalism derived in ref. |I| is correct, is 
was applied improperly in refs. |^, |||22l. Since our argu- 
ment rests on circumstantial evidence only, we strongly 
call for a remeasurment of the reaction ttN — > ujN close 
to the threshold. 

The paper is organized as follows. In Sec. Owe summa- 
rize the current status of the t:N-^ujN data and provide 
the data interpretation as it was given in the Ref. B . In 
Sec. [II we derive an analytical formula for the quantity 
measured in Refs. [H, 0, p[. In Sec. IV we describe the 
Monte Carlo simulations of the 7r~p— s-wn measurements. 
The evaluation of the imaginary part of the loN scatter- 
ing length from the irN^tuiN data is given in Sec. 0. 
The paper concludes with a summary. 



II. EXPERIMENTAL DATA 

In this section we summarize the current status of the 
data available for the ttN^luN total reaction cross sec- 
tion at pion beam energies close to the nominal pro- 
duction threshold, which is given by the pole mass 
M(^=781.94 MeV of the w-meson. There are three pub- 
lications from the same group that report results for 
T:~p—HLin measurements close to the nominal to produc- 
tion threshold |jl|, ||, ^ . The data of highest quality are 



presented in the paper of Karami et al. 0] , as shown in 
Table ||. The experimental procedure can be sketched as 
follows. 

The measurements were done at 33 different pion 
beam momenta p^ scanning the range from 1040 to 
1265 MeV/c. The beam momentum resolution was 
Ap^/pir— 0.8%, so each pion momentum setting addition- 
ally covered a range ~8-10 MeV. At each pion beam mo- 
mentum the neutrons were collected in 60 counters set 
within the angular range from 2.5° to 25.1° from the 
beam direction. 

The final statistics were regrouped in 2 MeV/c wide 
incident pion momentum bins with intensity about 10^ 
pions per bin. At each beam momentum the final neu- 
tron momentum (j„ and the neutron emission angle 9* 
in the center of mass system, as well as the fiux, were 
determined. The incident flux at all momenta was well 
fixed by weighting each event proportionally to the total 
number of the incident pions at its particular momentum. 

The data were distributed over 10 intervals of the neu- 
tron momenta q„ and 10 intervals of cosO* . This pro- 
cedure is equivalent to an integration over incoming tt- 
meson momenta for a fixed produced neutron momentum 
and angle. 

For every bin in qn and cos9* the missing mass spec- 
trum daexpl dm was reconstructed in order to separate 
the w-meson spectral distribution and the background 
|26| . The daexpl dm spectra were then fitted by a sum 
of an w-meson spectral function and a low-order polyno- 
mial. The fitting allowed separation of the background 
and w-meson signal. 

The spectral distribution da^xp/dm, at fixed intervals 
of g„ and cos0* was obtained by scanning or integrat- 
ing over the pion beam momentum p^. To replace 
the integration over p^ by an integration of the recon- 
structed missing mass spectra daexp/dm, over the mass 
771, the data were corrected by a corresponding Jacobian 
dp-^jdm. The Jacobian weakly depends on the beam mo- 
mentum. Neither in ref. ||3| nor in ref. ^ the inclusion 
of an additional Jacobian is mentioned. We come back 
to this point below. 

Finally the differential ir^p—^ajn cross sections were de- 
termined for all lOx 10 intervals in neutron momentum q„ 
and cos9* . The differential da/dfl cross sections shown 
in Ref. H for the different intervals of the neutron mo- 
menta are almost isotropic. The cross sections indicated 
in Table | were obtained by summing over the differential 
da/dfl cross sections for each neutron momentum inter- 
val. The qn intervals were specified for the central value 
of the cms neutron momenta P* given within the bin 
±AP/2. 

The data collected in that way were interpreted as 
two-body n^p-^ion reaction cross sections cr2b- These 
cross sections have been considered as those for a stable 
w-meson production, since the beam momentum inte- 
gration eliminates the dependence on the lu width. The 



TABLE I: The n^p^cun cross section a measured by Karaini 
at al. M for different intervals of the final neutron momenta 
P*±AP/2 in the center of mass system. Also are shown the 
reduced cross sections given by the ratio a/P* . 

J^ L\l (Texp ^exp / -t^ 

(MeV/c) (MeV/c) (fib) {pb/(MeY/c)) 



50 


20 


197±18 


3.94±0.36 


70 


20 


339±26 


4.84±0.37 


90 


20 


577±40 


6.41±0.44 


110 


20 


830±50 


7.55±0.45 


130 


20 


1118±71 


8.60±0.55 


150 


20 


1350±80 


9.00±0.53 


170 


20 


1510±74 


8.88±0.44 


190 


20 


1560±83 


8.21±0.44 



two-body reaction cross section is given explicitly by u4 
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where the factor of 6 accounts for the summation over the 
number of final spin states, the factor 2/3 is the isospin 
factor relating particle basis and isospin basis, M. is the 
spin averaged matrix element in the isospin basis and s 
stands for the squared invariant collision energy. Fur- 
thermore, in Eq. (^ (/„ and g„ are the incident and final 
center of mass momenta, respectively, with 



2^/i 
2^^ 



(4) 



where nin and jtitt are the nucleon and pion masses, re- 
spectively, Mi^=781.94 MeV is the pole mass of the w- 
meson and the function A is given by 



A(a;, y,z) = {x-y- zf - Ayz . 



(5) 



In line with Eq. (|^) , in order to extract the energy de- 
pendence of the matrix element squared jAlp one has 
to divide the measured cross sections by the phase space 
volume of the two-body final state. Already in the orig- 
inal publication the results on the reduced cross section 
a/P* (the last column of the Table ||) was given and the 
strong momentum dependence of the ir^n^nijp reduced 
cross section was interpreted as a strong momentum de- 
pendence of the reaction matrix element. 

The interpretation of the results given in Refs. ||l|, || ||] 
as a two-body reaction cross section was also used in the 
Landoldt-Bornstein compilation of elementary reactions 
cross sections H for both ir^p—^ujn and Tr~^n^ujp total 
cross sections. There the data are shown for fixed beam 
momentum corresponding to the center of mass momenta 
given in the experimental source. Note that the final 
center of mass momenta qn given in the experimental 
paper (here reproduced in Table m do not match the 



initial pion beam momenta given in the compilation |q| 
(here reproduced in Table ||) when calculated under the 
assumption, that a stable omega with its nominal mass 
is produced. How the values for p^^ in Rcf. [^ where 
determined is unclear to us. 



TABLE IL The n~p^Lun reaction cross section a and its 
uncertainty Act as given in the Landolt-Bornstein compila- 
tion M at the laboratory pion beam momenta p-^ range from 
1097 to 1170 MeV/c. Also shown are the matrix element 
squared \M\^ evaluated by Eq. (H) and the final neutron mo- 
mentum Qn in the center of mass system given by Eq. (Ml) . 



Ptt cross section error 
(GeV/c) CT (nb) Act ifib) 



\M\ 



qn 

(MeV/c) 



1.097 


197 


18 


182±17 


62 


1.102 


339 


26 


250±19 


78 


1.109 


577 


40 


347±24 


97 


1.118 


830 


50 


420±25 


116 


1.128 


1118 


71 


495±31 


135 


1.140 


1350 


80 


531±31 


154 


1.154 


1510 


74 


535±26 


174 


1.170 


1560 


83 


504±27 


193 



The spin averaged squared matrix element \M\'^ de- 
fined through Eq. (^) is shown in Fig.|l| as well as in Table 
y as a function of the neutron momentum g„ . The results 
are shown for the '!T~p—K^n (triangles) and n'^n—nup 
(squares) reactions. The data [g| on the 7r~p— s-wn re- 
action for laboratory pion beam momenta between 1.097 
and 1.17 GeV/c are shown by the circles. 

Fig.n^ illustrates the very strong momentum depen- 
dence of the squared matrix element near the reaction 
threshold at (7„<200 MeV/c mentioned above. Within 
the short range of the neutron momentum qn from 154 
to 62 MeV/c |A^p decreases almost by a factor of 3. 

Summarizing, there are three potential sources of en- 
ergy dependence of the TT~p^'UJn total cross section close 
to the nominal uj production threshold: the phase space, 
the width of the w-meson and the Tr~p—nun transition 
amphtude. It is the last quantity \M\'^ that contains all 
the tu production dynamics. In the energy regime of inter- 
est here the width effect on the energy dependence of the 
cross section is quite large. As was described above, to 
remove this effect from the data in a model-independent 
way, the authors of Ref.||l|, ||, |^ proposed to integrate 
over some range of initial pion momenta p^r for a fixed 
final neutron momentum labeled P*. It was shown in ref. 
||l| how the count rates can be related to the cross sec- 
tions. Due to the aforementioned difficulties for model 
studies to account for the resulting energy dependence 
of the cross section we repeat here the derivation from 
Binnie et al. and check the approximations made along 
the line with an event-by-event simulation. 
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FIG. 1: The spin averaged matrix element squared \M\ 
evaluated by Eq. (pi) from the total n~p—nun (triangles) 
and n'^n^'Ujp (squares) cross sections taken from the com- 
pilation M as a function of the final neutron momentum 
Qn in the center of mass system given by Eq. (m. The 
Tv~p~*Lon data quoted S] for the laboratory pion beam mo- 
menta 1.097<p,r<1.17 GeV/c are shown by the circles. 



III. ENERGY DEPENDENCE OF THE CROSS 
SECTION 

In Ref. Q it was argued that the energy dependence 
of the TT~p~*ujn reaction matrix element as shown in 
Fig. |l| can be traced to a misleading interpretation of 
the measured count rates. The basic findings were that 
it is indeed possible to extract the transition matrix el- 
ement independent of the shape of the mass distribu- 
tion of the resonance, however — in contradiction to the 
claims of Refs. [|| g, ^ — it was argued that the energy de- 
pendence of the count rates is different from that of the 
two-body reaction. As a consequence it was claimed, 
that the above mentioned energy dependence of the ma- 
trix element squared |A^p is a purely kinematical effect. 
In this section we want to check this claim by repeating 
the derivation of Binnie et al.. As it will turn out, the 
theoretical part of Ref. M is correct, however, there is a 
chance of misinterpretation of the final formula (formula 
(10) in Ref. O) which might have lead to a misuse of 
the technique. We will come back to this point below. 
In addition, we regard the following discussion as usefuU 
for future experiments for it is somewhat more compact 
than the original presentation of ref. [|l[ . 

Let us start with the cross section for ttN^XN reac- 
tion, where X denotes the decay products of the unstable 
meson. The reaction cross section can be expressed as 



da ^ {2-Kf 



9xA 



\M\HC, 



(6) 



As in Eq. (0) q.^ denotes the cms momentum of the initial 
pion and M. the spin averaged invariant matrix element. 
The factor of four is a combination of spin and isospin 
factors, as is given by Eq. (|3|). The trivial energy depen- 
dences of the reaction cross section are collected in the 
function C, which is defined as 

dC = d$/c+i(p;p„,Pi,-.,Pfc) \D{m^)W(j'i.-.Vk)\^ , 

(7) 
where we assume the unstable meson, whose propaga- 
tion is described by D(m}\ to decay into the fc par- 
ticles through the vertex function W . Here w? is the 
total invariant mass of the fc final decay particles with 
vr? = (X)-Pi)^j where the p^, i = \..k denote the corre- 
sponding four-momenta. Furthermore, the total initial 
four-momentum is denoted by p. The phase space of 
the final k particles and the final nuclcon with the four- 
momentum pn is defined as 

d^k+i{p\Pn,Pi. ■■■,Pk) = S'^{p-Pn ~ ^Pi) 

(27r)3 2^„li (27r)32w,' ^' 

where En and cui are the energy of the final nucleon and 
i-th decay particle, respectively. 

Using the unitarity relation we can introduce the spec- 
tral function p as 

(27r)3 fd<^kipx;Pi,-,Pk) \D{m'')W\'' = 



-QD{m ) =: p{m ). 



(9) 



A standard choice for the spectral function is that of a 
Breit-Wigner resonance. 



1 



p{m ) = — 



M^R 



^(™2_M2)2+M2r2 ' 



(10) 



where M^ and F;^ are the pole mass and width of the 
w-meson, respectively. It is this form that we will use 
in the evaluation of intermediate results. For the final 
result, however, the exact shape of the spectral function 
is irrelevant. All that we will use is the normalization 
condition, namely 



dm'^ p{m'^) = 1. (11) 

Now we can rewrite Eq. (m) in the center of mass system 



as 



^C = ./olxe Pi^"^) dm'^—^5{y/s - lo - E„) 
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1 
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■ p(s - 2^fsEn 



En 



(12) 



where u^^rn?- -1- q^^ is the energy of the final unstable 
meson state with the total invariant mass m and q„ de- 
notes the momentum of the final nucleon in the center of 
mass system. 



Finally, by introducing Eq. ( Jl^ ) into the Eq. (g), we 
obtain the differential reaction cross section as a function 
of the nucleon momentum: 



da 
dqn 



1 



27rg^- 



7i|l^l^^(- 



2^Eg,^+ml). (13) 



This formula agrees with that used in Ref. |g^ up to 
differences in the normalisation. Note that in the limit of 
vanishing width of the unstable particle Eq. ( |l3| ) trans- 
forms into the two-body cross section of Eq. p) for the 
production of particles with fixed masses, since 



lim p(s — 2\/sEn 



ml) 



EL 



2V^ U;. 



S{qn-q'„), (14) 



where E^ and q^ denote the energy and relative momen- 
tum of the produced stable particles. It is interesting to 
investigate for what values of F the right hand side of Eq. 
(|lj) is a good approximation to the left hand side. Nat- 
urally, the parameter that controls the behavior of the 
cross section should depend on the energy and the width 
of the decaying particle as well as the resolution of the 
detector. For the cross section to appear as a two-body 
cross section we should not resolve the decay particles 
any further. In Ref. El it was shown that for 



2P*AP 



>1 



(15) 



where /i denotes the reduced mass of the ujn system, the 
reaction rate behaves indeed like a two-body cross sec- 
tion. It should be obvious that the ratio AP/F appears 
here, for this ratio measures how closely we look at the 
production rates. For a given detector resolution AP we 
thus deduce a critical value for P* below which devia- 
tions from the two-body behaviour should be expected 
on purely kinematical grounds. Using AP — 20 MeV as 
given in Table | we find that, for values of P* > 90 MeV, 
the cross section for a fixed energy should behave like a 
two-body cross section. 

The solid lines in the Fig. H show the differential cross 
section da /dqn calculated for the T:^p—>ujn reaction at 
fixed TT-meson momenta p^r ■ Here we employed the rel- 
ativistic Breit-Wigner form as given in Eq. (HG) for the 
w-meson spectral function. Furthermore, in this calcu- 
lation, we employed the almost constant matrix element 
squared |A^P, as defined in Eq. (p4|). The neutron spec- 
tra are shown in Fig. 0for the (/„ range studied by Karami 
et al. g]. Fig. Q shows that the differential cross section 
da /dqn reflects the behavior of the spectral function. 

Fig. illustrates an important feature: namely, the 
differential cross section at different pion momenta, but 
at the same fixed neutron momentum are substantially 
different in absolute value. Naively one would expect 
that the da /dqn are identical for different p^ but the 
same qn- Note that the difference in the flux factor due 
to the different pion momenta is almost negligible within 
the range 1<Pt^<1.2 GeV/c. The calculations by Eq. (^) 
clearly show an absolutely different situation, since the 
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FIG. 2: The neutron momentum spectra calculated at differ- 
ent fixed pion beam momenta p,r. The histograms show the 
results from the event-by-even simulations, while the solid 
lines indicate the Eq. (M) . The spectra are shown only for the 
neutron momenta g„ range covered by Karami experiment |3| . 



phase space dependence of the cross section can not be 
factored from the spectral function density. 

It was the idea of Refs.|l|, ||, || to remove the depen- 
dence on the spectral function from Eq. (O) through an 
integration over the initial pion momentum p.,, . The nec- 
essary assumption is that the flux as well as the matrix 
element squared |A^p only weakly depend on the total 
collision energy when evaluated within a limited range of 
final qn momenta. In addition, we need to assume that 
neither of the two vary significantly when m is varied 
within the range where the spectral function is large. 

By integrating Eq. (g) over the laboratory 7r-meson 
momentum p^ we obtain the TT~p—nun cross section for 
the range of the neutron cms momentum P*±AP/2 as 



ct(P*,AP) 



dp-^ 



AP 



da 



1+ 

dqn , 
q_ dq. 



(16) 



where the limits of the integration over the neutron mo- 
mentum are fixed by 



q± = P*±AP/2, 



(17) 



while the integration over the pion laboratory momen- 
tum were performed 0] within the range from 1040 to 
1265 MeV/c. 
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FIG. 3: The differential cross section as a function of the pion 
beam momentum p^ calculated for different ranges of neutron 
momenta q„. The histograms are results from the event-by- 
event generator while the solid lines show the calculations 
using Eq. ([L6|). 

The solid lines in Fig. g show the neutron differential 
cross section as a function of the Tr-meson beam momen- 
tum calculated for the two different intervals of the neu- 
tron momenta. The shapes of the da/dp^^ distributions 
reflects the Breit-Wigner spectral function distribution. 
The width is dominantly given by the w-width with a 
tiny correction from the spread in the outgoing momen- 
tum. 

To remove the dependence on the spectral function p, 
we replace in Eq. (|1^) the integration over p^ by an inte- 
gration over the invariant mass squared m^ and addition- 
ally impose the normalization condition of Eq. (UW . As 
illustrated in Fig. the p^r integration can be assumed to 
be done over an infinite range. Therefore, the 7r~p— >wn 
cross section is given as 



a{P*,AP)^K 



with 



P*2 + -(APf 



dpn 



2T:Enqn\/s \drn? J 



\M\^ , (18) 



(19) 



where En=\ P* + vri^ and s denotes the invariant col- 
lision energy evaluated for the production of a stable w- 
meson with center of mass momentum P* . The weakly 



energy dependent Jacobian is given by 



yd-m? J AM^rUnqT, 



(20) 



where we have used the non relativistic expression for the 
total energy of the w-meson. 

The expression given in formula ( p^ ) agrees to eq. (10) 
of ref. |l|] for an ideal detector and beam, which states, 
translated into the quantities used above 



a = 2M^ 



pL)^ [ drd.^J^^ 



(21) 



where x^ denotes the cosine of the scattering angle in 
the lab frame, r denotes the time of flight of the outgo- 
ing neutron and Jr is the Jacobian relating the system 
(cos(6l),r) to (m,i) (c.f. eq. (6) in ref. 0) 



Jr 



{Ef^fM^d 



(22) 



where d denotes the flight distance for the neutrons and 
q^^{E^) denote the neutron momentum (energy) in the 
laboratory frame. In addition, for an isotropic matrix- 
element 



da 
'dt 



ttE,, 



^ \M\' 



dm'' 



dh 



dqnd^ 



IGnq'^s 



(23) 



Thus, at this stage it seems as if the experiments of refs. 
|ffl, 0, 0, Efl| where analysed properly. However, from 
what is written in refs. ||, |[ it is unclear, ii aexp, defined 
through 

a{P*,AP) = 2M^ (^\ ae.p{P*)- (24) 



is given, or da/dt as it can be extracted using eq. (|21 
The former option corresponds to interpreting Jt in eq. 
( ^ ) as J(T.q„) — the Jacobian that connects (cos(0),t) to 
{qn,t), where 



7 — J— T 

qn 



(25) 



where fl = M^Enj \fs^ which reduces to the reduced mass 
of the final state in the non-relativistic limit. This inter- 
pretation of the Jacobian appearing in eq. (pTl) is natural 
(though wrong) to the extend that the cut in the outgoing 
neutron momentum is given in terms of g„, whereas eq. 
( pl| ) is formulated in terms of the time of fiight r. The 
following evidence supports, that this misinterpretation 
acctually took place, namely 

• both publications [^, ^ only mention that there is a 
Jacobian to be included without further comment 
wether they used J-r (given in eq. ( P2|) or J(T,q„) 
(given in eq. (|25|), thus suggesting that only the 
naive one that connects the coordinate systems — 
namely J(r,g„) is used and 



in ref. 



there is explicitly only the analog of 



J(T,qn) mentioned and no other. Note, the momen- 
tum dependence of the amplitude given in ref. ^^ 
is consistent with that by Karami et al. [p7| . 

Thus one may consider two options: Either Jr (c.f. eq. 
(P2|)) was indeed included in the experimental analysis of 
Refs. ||l|, H, m ^ and the data. Then the strong energy 
dependence as shown in Fig. P is physical. 

Or the experimental data were analysed using J{T,q,^) 
as given by eq. (p5|). 

In what follows we investigate the consequences of the 
second option. In other words, we study the implications 
of the assumption that it was aexp defined in eq. (|24| ) 
that was given in the experimental papers [0, 0| . 

Matching Eqs. (^J) and (|l8|), we arrive at the central 
formula of this manuscript, providing the relation be- 
tween the measured cross section [|l|, g, ^ and the squared 
matrix element |A^p of the Tr^p—>ujn reaction, namely 



1 



^ exp 



AnqTrSp, 



P*' + {APy/12) \M\' 



(26) 



One might wonder why the resolution of the neutron de- 
tector enters the final expression. On a second thought 
this should be quite obvious, since for every given value of 
AP there will be a count rate even for a minimal value of 
P* . This is most clearly demonstrated in the uppermost 
panel of Fig. @). 

By construction, A4 is the matrix element for the reac- 
tion TT~p^ojn considering the w-meson as a stable par- 
ticle, since the width of the omega was treated explicitly. 
We can therefore define an effective two-body cross sec- 
tion a'^-^-f given by Eq. (p|) in terms of the cross section 
aexp measured in Ref. [ 



yeff 



fiP* 



P*' + {AP)^/12 



(27) 



Therefore, the effective cross section a'^^^ is the 
Tr~p—HL!n total reaction cross section taken under the 
assumption that the produced a;-meson is s table with 
the fixed mass M^=781.94 MeV. In Table |n| we provide 
the corrected cross section a'^^^ based on the hypothesis 
that the formalism given in Ref. [y was applied im- 
properly. The corrected Tr~p—fUJn reaction cross sections 
deviate from the results quoted in the compilation [g| by 
a factor up to 8 within the range of pion beam momenta 
1094<p^<1167MeV/c. 

We can now evaluate the corrected matrix ele- 
ment squared \Ai\'^ from the Tr~p—>u!n data employing 
Eq. (pi). The results are shown in Fig. and clearly il- 
lustrate that the corrected matrix element is in line with 
the other data and does not indicate any pathological 
energy dependence at small neutron momenta. The re- 
sulting energy dependence is in a lot better agreement 
with that of model studies pS] . This could be taken as 
further evidence, that indeed the formalism derived in 
Ref. m was applied improperly. 



TABLE III: The corrected effective Tr'p^ujn cross section 
a'^-'-' as derived by us. The second to last column contains the 
error as it is derived directly from the data H , whereas in the 
error given in the last column the uncertainty of the neutron 
momentum P* is included as well. The relation between the 
central neutron momentum P* and the pion beam momentum 
p,r is given by Eq. (W). 



P* p, a"" 
(MeV/c) (MeV/c) (pb) 


Sagf 
(/.fe) 


{^^b) 


50 1094 


1656 


151 


364 


70 1099 


2049 


157 


332 


90 1106 


2720 


189 


356 


110 1115 


3206 


193 


350 


130 1125 


3656 


232 


365 


150 1137 


3828 


227 


341 


170 1151 


3780 


185 


289 


190 1167 


3494 


186 


262 



IV. EVENT-BY-EVENT SIMULATIONS 

In addition, to check the analytical calculations of the 
previous section, we also developed an event-by-event 
generator to simulate the experimental measurements de- 
scribed in Refs. 1^, g, ||. In addition, this simulation 
allows us to discuss several aspects of the technique de- 
veloped in Ref. [Q in detail. 

In accordance with Ref. H, the Tr~ p^wn events were 
simulated at pion momenta randomly selected within the 
range between 1040 and 1265 MeV/c. At fixed pion mo- 
mentum Ptt the events were generated by the following 
method. 

1. The flux factor was determined as 



C = 8Tr^X^^^{s,ml,ml), 



(28) 



where rrijr, rup and to„ are the pion, proton and neutron 
masses, respectively and s is the squared invariant energy 
of the event given by 



= rn^+ml + 2mn\J pi ^ m\. 



(29) 



2. The squared oj-meson mass rr? was randomly gen- 
erated by the Breit-Wigner probability distribution as 
defined in Eq. (|l^). 

3. The neutron momentum was calculated as 



A-^/^(s,m^,TO^J 

and the phase space factor of the event is given as 
A-^/^(s,m^,TO^) 



$ = TT 



2s 



(30) 



(31) 



4. Now each event at fixed pion beam momentum p-^ 
and with certain w-meson mass to was accounted for with 
the weight (pi 



VFe,„ 



^\Mf 



(32) 
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FIG. 4: The spin averaged matrix element squared \M\'^ 
extracted from the total 7r~p— >ajn (triangles) and Tv^n—nop 
(squares) cross sections taken from the compilation |5| as a 
function of the final neutron momentum q„. The circles show 
the corrected matrix element squared evaluated using Eq. (Eq) 
for the data |3| on the Tv~p^Lon reaction. The solid line shows 
the parametrization of the matrix element given by Eq. ( M) 
with the parameters fixed by a fit of the event-by-event sim- 
ulation of the data of Ref. [HI . 



As in Eq. (g|), Ai denotes the spin averaged matrix el- 
ement in isospin basis and the factor of 4 stems from a 
combination of the number of final spin states and the 
isospin factor. Note that this expression agrees with that 
for a two-body cross section (c.f. Eq. (||)), however, in 
our analysis we treat the energy and mass dependence of 
Qn properly. 

5. In order to ensure energy conservation we impose 
the kinematical condition Wev=0 when ^/s<m+mn. 

Now, to calculate the 7r~p— s-wn cross section at fixed 
pion momentum a{pTr), we generate N events and sum 
them as 



ev=N 



<p-) = ^ E w^^ 



(33) 



As described in the previous section, in Ref. P] the 
data are given not at fixed pion beam momenta, but at 
different ranges of the neutron momenta. At fixed pion 
beam momentum the reaction TT~p—nvn produces the 
neutron momentum spectrum rather than a fixed neu- 
tron momentum because of the variation of the w-meson 
mass. It can be well understood within the event-by- 
event simulations when proceeding from step 2 to step 3, 
as illustrated in Fig. |[ The histograms in Fig. || show 
the results from the Monte Carlo simulations, which are 
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FIG. 5: The missing mass spectra calculated for the fixed 
pion momenta p^. The full histograms show the results 
from the event-by-even generator. The solid lines indicate 
the Breit-Wigner distribution for the cj-meson. The hatched 
histograms show the mass spectra for the neutron momenta 
ranges 60<q„<80 MeV/c (A) and 140<<7„<160 MeV/c (B). 
The excess energy e = y/s — nin — M^j . 



in agreement with the analytical evaluations. 

Fig. also shows the missing mass distribution sim- 
ulated for the Tr~p—fLon reaction at different fixed pion 
beam momenta. The histograms in Fig. |5| indicate the to- 
tal missing mass spectra, while the solid lines indicate the 
relativistic Breit-Wigner distribution for the a;-meson. 
It is clear that at low pion momenta only w-mesons with 
small masses can be produced because of the energy con- 
servation imposed by the condition of step 5. 

The hatched histograms in Fig. ||, moreover, show 
the missing mass spectra for the two neutron momenta 
ranges 60<g„<80 MeV/c (A) and 140<g„<160 MeV/c 
(B). Therefore, for a fixed range of final neutron mo- 
menta, the contribution from the different w-meson 
masses comes due to the different incident pion momenta. 
Fig. most clearly illustrates the basic idea of the mea- 
surements [y, y, 1^ : namely to saturate the cj-meson spec- 
tral function by scanning the pion beam momentum. As 
indicated by the hatched histograms, the lu production 
cross sections at different 7r-meson momenta but at the 
same neutron momentum qn substantially differ in abso- 
lute value. 

An additional generation over the inital pion momen- 
tum Ptt is required, randomly scanning the experimental 



pion momentum range. After introducing the Jacobian 
defined in Eq. (po|), we fit the experimental data [|| in or- 
der to extract the matrix element squared |A^p as well as 
its energy dependence. For the squared matrix element 
|A^P we assume an energy dependence given by 



1-^1' = TX^' 



(34) 



with parameters Mq and b to be fit to the data of Ref. 

Fig. g shows the experimental results from Ref. to- 
gether with the results of the fit with the parameters 

M^ = (2.1 ± 0.1) X 10^ and b = (0.7 ± 0.1) fm^ . (35) 

The experimental cross sections aexp pl are shown by the 
solid circles, while the boxes are our calculations. The 
size of the boxes shows the uncertainty in the extraction 
of the matrix element squared \M\'^ as well as the range 
of integration used for the neutron momentum as given 
by AP. 

Finally, the squared matrix element parametrized by 
Eqs. (p4) and (p3) fitted by using the Monte Carlo simu- 
lations to the experimental data [gj is shown by the solid 
line in the Fig. H. It agrees well with the results derived 
analytically. We also find that the parametrization given 
by Eq. ( p4| ) works rather well over a wide energy range. 
Note that the fit was carried out on the basis the dataset 
of Ref. H only, thus only the momentum range from 50 
to 200 MeV was fit. 



V. EVALUATION OF THE SCATTERING 
LENGTH 

The matrix element at vanishing neutron momentum 
is a quantity of particular interest. In order to evalu- 
ate I Alp at qn=0 we fit the results of the Fig. [^ with a 
polynomial in qf^ as suggested in Ref. 0| 



\M\^=a + bql 



cqn 



We find for the corrected matrix element a = 2.0 
b = —0.6 fm^ and c = —0.3 fm"* and therefore 



lim \MY 



(2.0 ±0.4) X 10^ 



(36) 



10^ 



(37) 



It is reasuring that the different parametrization of Eq. 
( p4| ) lead to the same value of the matrix element at 
threshold within the experimental accuracy. 

It is now straight forward to extract a lower limit for 
the imaginary part of the elastic scattering length. As- 
suming that the ujN—^ttN reaction channel provides the 
most significant part of the inelasticity for the coN scat- 
tering, the value derived should be a reasonable estimate 
of the true imaginary part of the scattering length. Note 
that model calculations O, nj] show that the nN chan- 
nel indeed dominates the ujN interaction. 
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FIG. 6: The cross sections a^xp measured at Karami experi- 
ment |m] at different intervals of a neutron momenta q„. The 
experimental results M are shown by solid circles, while the 
boxes indicate our calculations. The size of the boxes illus- 
trates the uncertainties in evaluation of the reaction matrix 
element \A4\ and the neutron momentum interval. 



We start with the S'-matrix describing meson-baryon 
scattering near the w-meson production threshold and 
assume n physical final state channels to be relevant at 
this energy. The presentation is given in the isospin ba- 
sis. Let the ujN elastic channel be first, while the irN is 
second. The matrix element corresponding to the elastic 
LuN channel is parametrized as 



5ii :— r/exp(2i(S), 



(38) 



with r] being the inelasticity in the ujN scattering and 6 
is the elastic scattering phase shift. 

The unitarity constraint SS'' = 1 now translates into 
the relation 



i=2 



l)|5lJ 



(39) 



where the factor {2si + 1) contains the spin multiplicity 
of the particular channel i. Since all the contributions in 
the sum arc positive, Eq. (39) directly leads to 



l-?7^>2|S'i 



(40) 



For s-wave scattering the relation between the 5*- 
matrix and A4 evaluated previously is given by 



XiXj 



2(27r)3 



ujiEiUjjEj 



(47r)Al 



i-j, 



(41) 
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where 



TT —Pi 



V~s 



(42) 



with Pi being the on-shell momentum in the cms of the 
particular system and Ei and coi denoting the baryon and 
meson energy in that channel, respectively. Combining 
the two last equations leads to 



l-V > 



^ qn^\M\', 



(Att) 



(43) 



where 5^ and (/„ are the incident pion and final nucleon 
momenta in the center of mass system, respectively. 

The scatting amplitude f^N in the ujN channel is given 
by 



fu 



■N 



2qni 



1 



776 



2iS\ 



and thus its imaginary part reads 



1 



^Un = ~-z — (1 - ?7cos(2(5)). 



(44) 



(45) 



The low momentum behavior of both the inelasticity 
{1 — rj) as well as the phase shift 6 is linear in the mo- 
mentum. Therefore, close to the nominal uj production 
threshold we evaluate the imaginary part of the coN scat- 
tering length SfltjAr as 



1 



^OwAT 



lim 3/^jv > „., ^2 



Qn 



\M\' 



(46) 



Note that this formula agrees with the one given in 
Ref. 0. 

Finally we deduce for the imaginary part of the spin 
averanged scattering length 



3a^7V > (0.24 ±0.05) fm. 



(47) 



As can be seen from Fig. M, the value for the corrected 
squared matrix element at threshold changed by more 
than an order of magnitude. Since the lower limit for 
'^o.ljN scales linearly with the squared matrix element, 
our bound naturally is significantly stronger than those 
present in the current literature. For instance, in Ref. |1^ 
a value 9acjAr>0.02 fm is given (Note: in Ref. Q the 
scattering length was deduced including the iT~p channel 
only. To include the tt^tt, isotopic channel the result of 
Ref. ua was multiplied with the isospin factor 3/2 to 
allow comparsion to our value) . 
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FIG. 7; The matrix element extracted from the data of 
Ref. ly. The opaque circles denote the corrected matrix ele- 
ment using Eq. (ETI) whereas the filled circles show the uncor- 
rected result. The small error bars are those from the data 
only, whereas the large error bars include the uncertainty in 
P* , given by AP=20 MeV. The solid lines are polynomial fits 
to the data. 



improperly applied in the analyses. A purely kinemat- 
ical factor not only corrects this error and removes the 
unusual energy dependence resulting from the primary 
analysis, but also keeps the data in accordance with the 
world data set. 

Based on the modified data we extracted a value for a 
lower bound the imaginary part of the elastic coN scat- 
tering length. The new value is larger by an order of 
magnitude compared than that used in the literature so 
far ||ll|. As should be clear from the discussion in the 
introduction, a large change in the imaginary part of the 
scattering length should have a large effect on estimates 
for both the in medium width of the omega as well as the 
existence of w-nucleus bound states. 

We would like to stress that we call for a remeasure- 
ment of the ttA^ —^ ldN reaction in the kinematics close 
to the nominal to threshold given the importance of this 
reaction, as outlined in the introduction. This remea- 
surement will not only allow to confirm the claims of 
this work but should also fill the gap in the experimental 
data between 200 and 500 MeV/c final center of mass 
momenta. 
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VI. SUMMARY 

Using a Monte Carlo simulation as well as an analyt- 
ical calculation we reinvestigated the experimental data 
of Refs. |0, pi H] Based on circumstantial evidence, we ar- 
gued, that the formalism derived correctly in Ref. ||l| was 
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